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We  present  numerical  simulations  of  \elocity  modulated  held  effect  transistors  as  proposed  bv  Sakaki.  Using  self-eon- 
sistent  particle-held  Monte  Carlo  analysis,  we  assess  possible  advantages  of  these  novel  device  structures  with  respect  to 


switching  speed  and  show  the  qualitative  correctness  ol  Sa! 
redistribution  of  electrons  in  the  separated  channels. 

§1.  Introduction 

In  field  effect  transistors,  the  source  to  drain  current  is 
controled  by  modulating  ihe  charge  eauier  concentration 
under  the  gate.  Therefore,  the  ultimate  switch-on  time  is 
roughly  determined  by  the  transit  time  of  electrons  from 
source  to  drain.  The  calculated  switch-on  time  for  a  0.4 
micrometer  gate  length  (I  micrometer  from  source  to 
drain)  high-electron  mobility  transistor  is  longer  than  3.5 
picoseconds  at  300  K,i:'  and  about  3.2  picoseconds  at  77 
k. 

To  achieve  sub-picosecond  switching  in  semiconduc¬ 
tors,  Sakaki  proposed  the  concept  of  velocity-modula¬ 
tion  transistois  (VMI’s)."  The  velocity-modulation 
transistor  concept  attempts  to  capitalize  on  the  extremely 
short  perpendicular  transit  times  between  two  adjacent 
channels,  with  different  transport  properties  (i.e,,  mobili¬ 
ty).  The  original  structure  suggested  by  Sakaki  is  a  n- 
AIGaAs/GaAs/n-AlGaAs  double  heterostructure  (Fig. 
U."  Two  channels  are  formed  at  the  heterointerfaces  in 
GaAs.  One  of  the  channel*  is  heavily  doped  w  ith  accep¬ 
tors  and  donors,  and  therefore  exhibits  a  low  mobility, 
while  the  other  channel,  the  high  mobility  channel,  is  un¬ 
doped.  Schottkv  barriers  are  placed  on  top  and  bottom 
of  the  device  to  switch  the  electrons  back  and  forth 
perpendicular  to  the  interface,  hence  modulating  the 
mobility  and  current  between  the  source  and  drain. 

If  it  is  assumed  that  the  intrinsic  device  speed  is  limited 
by  the  perpendicular  transit  time,  the  switching  time  can 
be  estimated  as  follows.  For  a  channel  separation  of  500 
A,  and  a  velocity  of  1  x  10  cm/s  (ignoring  velocity)  over¬ 
shoot)  perpendicular  to  the  heterointerface,  the  swit¬ 
ching  time  can  be  estimated  as  r  =  0.5  picoseconds.  This 
is  an  order  of  magnitude  faster  than  common  switching 
times  which  involve  the  gate  length.'4’ 
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l  ig.  1.  Device  topology  of  the  velocity -modulation  transistor  (YVT) 
as  proposed  hv  Sakaki. 


aki‘s  ideas.  Quantitatively,  delays  are  introduced  by  the 


The  dynamics  of  electron  transfer  between  two  adja¬ 
cent  equivalent  high  mobility  channels  (dual-channel 
high-electron  mobility  transistor)  were  previously 
reported  for  300  k  operation.^1  It  has  been  concluded 
that  the  ideas  of  multiple  channel  transfer  are  valid,  and 
ultimate  device  speeds  could  be  enhanced  if  these  ideas 
can  be  realized.  However,  it  also  has  been  found  that  the 
electron  redistribution  along  the  device  channel  plays  an 
important  role  and  prolongs  the  time  period  which  is 
necessary  to  reach  steady  stale,  even  if  the  electrons  are 
switched  perpendicular  to  the  layers  rather  fast. 

The  purpose  of  this  paper  is  io  discuss  the  factors  in¬ 
fluencing  the  t i uK  constants  of  a  velocity-modulation 
transistor,  at  77  k  and  300  k.  We  will  elaborate  on  the 
ultimate  device  speed,  on-state  and  off-state  current 
levels,  time  evolution  of  the  electron  distribution  func¬ 
tion,  and  ionized  impurity  scattering  as  a  means  to 
modulate  channel  mobility  of  a  field  effect  transistor. 

§2.  Device  Structure  and  Simulation 

The  structure  used  in  the  simulations  in  an  n-AIGaAs/ 
GaAs/n-AlGaAs  double  heterostrueture  as  shown  in  Fig. 
2.  The  AlGaAs  layers  are  doped  relatively  lightly 
(/Vj)  =  2  x  10'  cm  ',  where  is  the  donor  concentration). 
The  two  channels  (A  and  B)  formed  at  the  heterointer¬ 
face  are  assumed  to  be  contacted  separately.  Channel  B  is 
assumed  to  be  undoped  (,\j,  =  1  x  1 0'‘  cm  ')  and  channel 
A  is  heavily  doped  with  both  acceptors  and  donors  pre¬ 
sent  in  order  to  reduce  the  mobility.  Schottkv  barriers  are 
placed  on  lop  and  bottom  of  the  device  to  modulate  the 
conductivity  ol  the  channels.  The  thickness  of  the  GaAs 
is  chosen  to  be  1000  A.  Under  these  circumstances  size 
quantization  effects  are  unimportant  and  our 
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semidassieal  treatment  is  appropriate.  The  source  to 
drain  distance  is  1  micrometer  and  the  gates  are  sym¬ 
metrically  placed  between  source  and  drain  and  are  0.4 
micrometers  long. 

Our  Monte  Carlo  model  incorporates  the  T-L-X  band 
structure  for  both  GaAs  and  AIGaAs.  Polar  optical 
phonon  sattering,  equivalent  and  non-equivalent  inter¬ 
valley  scattering,  real  space  transfer,  impact  ionization 
and  ionized  impurity  scattering  are  included  in  the  model 
which  has  been  described  in  detail  previously.^"' 

The  total  number  of  sample  electrons  in  the  device  is  a 
variable,  and  the  electron  population  is  allowed  to  adjust 
itself  to  satisfy  overall  charge  neutrality.  Small  regions 
surrounding  the  ohmic  contacts  are  kept  neutral  by  injec¬ 
ting  electrons  chosen  from  an  equilibrium  energy  distribu¬ 
tion  as  needed.  The  source  and  drain  currents  are  deter¬ 
mined  by  recording  the  number  of  injected  and  absorbed 
particles  during  a  preset  time  interval.  M 

The  electric  field  distribution  in  the  device  domain  is 
updated  every  time  step  (2  femtoseconds)  by  solving 
Poisson’s  equation  numerically  by  the  finite  difference 
method  and  successive  over  relaxation.  Dirichlet  condi¬ 
tions  are  imposed  at  the  metal-semiconductor  surfaces, 
and  potential  continuity  conditions  are  used  at  the 
AIGaAs/GaAs  interface.  The  perpendicular  field 
changes  as  the  ratio  of  the  dielectric  constants. 

§3.  Results 

3.1  Switching  times 

We  first  present  the  switching  time  simulations  for  300 
k  operation.  The  dev  ice  is  initially  (t  =  0  )  in  its  oil-state 
as  shown  in  Fig.  3(a).  The  current  is  carried  by  electrons 
in  the  low  mobility  channel  (channel  A).  The  gate  voltage 
of  channel  A  is  0.6  volts,  and  the  gate  voltage  of  channel 
B  is  0.0  volts.  Source  voltages  are  kept  at  0.0  volts  and 
the  drain  voltages  are  kept  0.8  in  all  of  the  simulations. 
To  switch  the  device  on,  the  gate  voltage  of  channel  A  is 
instantaneously  decreased  to  0.0  volts,  and  the  gate 
voltage  of  channel  B  is  instantaneously  increased  to  0.6 
volts.  Channel  A  is  immediately  fully  depleted  (Figs.  3(b) 
and  3(c)).  The  transit  time  for  perpendicular  transport 
(from  one  interface  to  the  other)  is  completed  in  about 
0.2  picoseconds,  regardless  of  the  doping  profiles. 
However,  steady  state  is  only  completely  reached  after 
the  electrons  that  are  transferred  to  the  high  mobility 
channel  (channel  B)  redistribute  themselves  to  achieve 
normal  channel  conduction.  This  redistribution  takes  2.0 
to  3.0  picoseconds,  depending  on  the  ionized  impurity 
_  concentration  of  the  low  mobility  channel  ((Fig.  3(d)). 

The  switch-off  simulation  starts  from  the  steady  state 
on-configuration  (Fig.  3(d)).  The  gate  voltages  of  chan¬ 
nel  A  and  B  are  instantaneously  reversed  as  shown  in 
Figs.  4(a)  and  4(b).  The  electrons  are  transferred  to  chan¬ 
nel  A  again  in  about  0.2  picoseconds.  Howevey  steady 
state  is  reached  only  after  2.5  to  2.8  picoseconds  depend¬ 
ing  on  the  ionized  impurity  concentration  in  the  low 
mobility  channel  (Figs.  5(a)— 5(d)),  again  because  of  the 
redistribution  of  electrons  in  the  given  channel. 

We  have  defined  steady  state  as  the  point  of  operation 
where  the  drain  and  source  currents  have  stabilized  and 
are  within  100<j  of  each  other.  After  0.2  picoseconds. 
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I  ig.  .V  hlectron  transfer  behavior  during  switch-on  transient  at, 
(a)  /-()  picoseconds,  (b)  /=■().  I  picoseconds,  (c)  /“0.2  picoseconds, 
(d)  /  =  3  picoseconds. 


which  is  the  perpendicular  transit  time,  source  and  drain 
current  levels  can  be  a  factor  of  two  apart  from  each 
other,  or  their  respective  steady  state  values.  We  con¬ 
clude  therefore  that  the  intrinsic  device  speed  is  not  only 
limited  by  the  perpendicular  transit  time. 

The  switch-on  and  switch-off  times  w  hich  are  necessary 
to  reach  steady  state  at  3C0  K  and  77  K  for  various  dop¬ 
ings  of  channel  A  are  summarized  in  Table  I.  The  switch- 
on  and  switch-off  times  calculated  for  the  velocity- 
modulation  transistor  (VM'ij  are  siiortei  than  for 
longitudinal  (source-to-drain)  switch-on  times  which  are 
3.5'  and  3.2  picoseconds  for  300  k  and  77  k  respec¬ 
tively. 

Increasing  the  impurity  concentrations  leads  to  longer 
switch-on  and  switch-off  times.  The  calculations  also 
show  that  switch-off  is  always  faster  than  switch-on.  We 
comment  on  this  point  in  the  discussion  section. 


Ensembk  A U 


Carlo  Simulation  of  a  Velocity- Modulation  Held  Effect  Transistor  (l  \  IT) 


1521 


t=0.1ps 

Gate  A 


Gate  B 


lb) 


t  =  0  2  ps 


Gate  A 


Gate  B 


(d) 

tig.  5.  Electron  transfer  behavior  during  switch-oil'  transient  at, 
la)  r  =  0  picoseconds,  (b)  r  =  (),l  picoseconds,  (cl  /=(). 2  picoseconds, 
(d)  (  =  3  picoseconds. 


Table  1. 

Channel  A  doping  r„„  (300  K)  r.„,  (77  K)  r,„,  (3(H)  K)  r,M  (7"1  K) 


Tig.  4.  The  energy  band  profile  of  the  structure  across  the  gates 
(at  x=0.5  micrometers)  for  two  different  biasing  conditions, 
(a)  T(lA  =  0.0  V,  T„„  =  0.6  V,  (b)  )„v=0.6V,  T(,„  =  0.0V. 

3.2  Evolution  of  the  electron  distribution  function 
The  election  energy  distribution  function  cannot  di¬ 
rectly  be  measured  but  is  important  for  the  under¬ 
standing  of  the  device.  The  study  of  its  time  evolution 
shows  clearly  that  the  intrinsic  device  speed  is  not  only 
limited  by  the  perpendicular  transit  time  (0.2  pico¬ 
seconds),  but  also  by  the  lateral  distribution  of  electrons. 

For  the  biasing  conditions  considered  in  this  paper,  the 
electrons  in  the  device  channel  are  mostly  in  the  T -valley, 
except  at  the  drain  edge  of  the  gale  (pinch-off  region). 
Therefore  only  the  T-vailey  electrons  in  GaAs,  under  the 
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gate  region  of  the  device  will  be  described  in  detail.  - 

Initially  (/  =  0“ ),  the  electrons  are  in  the  low  mobility 
channel,  and  are  at  steady  state  (Fig.  6(a)).  The  displace¬ 
ment  of  the  distribution  function  in  the  A\  direction  - 

(parallel  to  the  interface)  indicates  some  current  flow 
from  source  to  drain  in  the  device.  When  the  gate  biases  or 
are  reversed  ( F<,s  =  0.0jajit,  F,.R=0.6  volts  at  /  =  ()' ),  the 
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Fig.  6.  Evolution  of  the  electron  distribution  function  under  the  gates  during  switch-on  transient  at,  (a)  /=■  0  picoseconds, 
(b)  /  =  0. 1  picoseconds,  (c)  r=0.2  picoseconds,  (d)  /  =  3  picoseconds. 


electron  distribution  function  is  shifted  into  the  —k, 
(perpendicular  to  the  heterointerface)  direction 
(Fig.  6(b)).  Electrons  impinging  on  the  interface  are 
reflected  back  into  GaAs  if  they  cannot  surmount  the 
potential  barrier.  Therefore  the  first  quadrant  of  A-space 
is  populated  (Fig.  6(c)).  Steady  state  is  reached  sometime 
after  the  electrons  have  transferred  to  the  high  mobility 
channel  (Fig.  6(d)). 

The  evolution  of  the  electron  distribution  function  dur¬ 
ing  the  switch-off  transient  is  shown  in  Figs.  7(a)-7(d). 

It  is  observed  that  the  transient  and  steady  state 
distribution  functions  are  not  spherically  symmetric. 
This  deviation  from  a  Maxwell-Boltzmann  distribution 
follows  from  the  dependence  of  the  scattering  rates  on 
the  particle  momentum,  and  energy  exchange  between 
electrons  and  polar  phonons  in  GaAs. 

3.3  Current  levels 

In  the  operation  of  the  velocity-modulation  transistor, 
an  increase  in  source-to-drain  current  is  expected  as  the 
electrons  are  moved  to  the  high  mobility  channel.  The 
current  switch-on  transient  at  77  K  is  shown  in  Fig.  8. 
The  off-state  current  is  /<»i  =32  mA/mm  (at  /  =  0).  After 
0.2  picoseconds  (perpendicular  transit  time  of  electrons) 
the  source  and  drain  currents  have  not  yet  reached  steady 


state,  which  is  reached  in  =2.4  picoseconds  and  the  on- 
state  current  is  /,)N  =  98  mA/mm. 

At  300  K,  a  current  of  /()N  =  80  mA/mm  flows  through 
the  VMT  when  it  is  in  on-state.  The  “off-state"  current 
is,  /on  =49  mA/mm  for  a  heavily  doped  low  mobility 
channel  (A\  =  2x  10'' cm1,  /V„  =  2.001  x  10'"  cm  '). 

Therefore,  the  channel  current  can  be  modulated  by 
=  40%  at  300  K,  and  by  =70%  at  77  K,  if  a  highly  com¬ 
pensated  concentration  of  ionized  impurities  is  introduc¬ 
ed  into  the  device  channel. 

3.4  Average  electron  velocity  distribution 

The  average  electron  velocity  along  the  device  channel 
is  sketched  for  300  K  in  Fig.  9,  and  for  77  K  in  Fig.  10  for 
various  channel  dopings. 

The  velocity  distribution  profiles  for  all  channel  dop¬ 
ing  have  the  same  qualitative  form.  The  electrons  enter¬ 
ing  the  high  electric  field  region  at  the  drain  side  of  the 
gate,  experience  velocity  overshoot.1"  The  peak  electron 
velocities  are  4.1  x  10’  cm/s  at  300  K,  and  6.3  x  10  cm/s 
at  77  K.  The  average  electron  velocities  in  the  low  electric 
field  regions  of  the  device  are  much  lower.  The  effect  of 
the  ionized  impurities  on  the  velocity  distribution  profile 
is  more  prominent  at  lower  temperatures,  especially  in 
the  low  electric  field  regions  (i.e.,  source-to-gatc)  of  the 
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Fig.  7.  Evolution  of  the  electron  distribution  function  under  the  gates  during  switch-ofT  transient  at.  (a)  /=()  picoseconds, 
(b)  r  — 0.1  picoseconds,  (ct  t—  0.2  picoseconds.  <d )  /  =  3  picoseconds. 


Fig.  8.  Current  switch-on  transient  for  the  velocity-modulation  tran¬ 
sistor  at  77  K. 


Fig.  9.  Average  electron  velocity  in  undoped  and  heavily  doped  chan¬ 
nels  of  the  velocity-modulation  transistor  at  300  K. 
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l  iji.  10.  Average  electron  velocity  in  uiuloped  and  hcav  ily  doped  chan¬ 
nel*  of  the  velocity-modulation  transistor  at  77  k. 


device.  This  is  due  to  the  increased  ioni/ed  impurity  scat¬ 
tering  rate  at  lower  temperature  and  lower  average  elec¬ 
tron  energies. 

§4.  Discussion 

The  calculations  presented  above  show  that  the  current 
through  a  device  can  be  reduced  by  40°o  at  300  k  and 
70%  at  77  k  if  electrons  are  transferred  from  a  high 
mobility  channel  to  a  highly  doped  compensated  channel 
with  <V„  =  2.001  x  10l!i  cm  '  and  ,V\  =  2x|0''cm  ’. 
However,  the  reduction  in  electron  velocity  distribution 
in  the  device  under  the  gate  is  at  most  30%  at  300  k,  and 
60%  at  77  k.  Therefore  velocity  reduction  alone  does  not 
explain  the  total  change  in  channel  conductance.  An  in¬ 
vestigation  of  the  electron  distribution  in  the  conducting 
channels  of  the  device  (Fig.  3(a)  and  Fig.  5(a))  reveals 
that  there  are  more  carriers  in  the  channel  when  the 
device  is  in  its  on-state  (Fig.  5(a)).  In  fact,  the  concentra¬ 
tion  of  electrons  in  the  channel  is  increased  by  about 
15%  when  the  device  is  switched  on. 

Change  in  the  carrier  concentration  under  the  gate  is 
brought  about  by  the  different  source  resistances  and 
field  distributions,  depending  upon  whether  the  electrons 
are  in  undoped  or  heavily  doped  channel.  The  difference 
in  the  source  resistances  can  be  explained  from  the  energy 
dependence  of  impurity  scattering.  As  it  is  well  known, 
ionized  impurity  scattering  prefers  small  scattering 
angles,  especially  at  higher  electron  energies.1"1  Therefore 
ionized  impurity  scattering  is  most  important  in  determin¬ 
ing  mobilities  in  the  low  electric  field  regions  of  the 
device,  i.e.,  the  source  to  gate  region.  This  difference  in 
source  resistances  is  sufficient  to  explain  the  higher  elec¬ 
tron  concentrations  under  the  gate  when  the  channel  is 
undoped. 

Finally  we  would  like  to  comment  on  the  observed 
asymmetry  between  switching  on  and  off.  The  calculated 
switch-on  times  are  always  larger  than  switch-off  times. 
This  is  a  general  property  of  the  field  effect  transistors 
and  in  this  case  also  a  consequence  of  the  symmetric  gate 
biases  used  for  switching  here.  As  mentioned  above, 
change  in  the  channel  conductance  is  due  to  two  factors; 
change  in  electron  mobility  and  electron  concentration 
{AG  =  nA fi  +  fiAn,  where  AG  is  the  change  in  channel 


conductance,  J  g  is  the  difference  in  mobility  between  the 
channels,  and  An  is  the  difference  in  electron  concentra¬ 
tion  under  the  gate).  To  switch  the  VMT  on,  by  simply 
reversing  the  gate  biases  will  bring  in  more  electrons 
( —  15 ®o  more  than  what  has  been  transferred  from  the 
low  mobility  channel)  from  the  source  and  drain  contact. 
To  switch  the  VMT  off,  electrons  are  transferred  to  the 
low  mobility  channel  and  any  excess  electrons  are  dispos¬ 
ed  of  through  the  source  and  drain,  i.e.,  excess  carriers 
need  not  transit  the  total  source  to  drain  (or  at  least  gate) 
length. 

This  suggests  that  switch-on  and  switch-off  times  can 
be  reduced  further  (to  about  2  ps  at  300  k)  if  slightly 
asymmetric  gate  biases  which  force  the  condition  J/t  =  0. 
for  the  operation  of  the  VMT  are  used.  However,  this 
will  be  at  the  expense  of  a  smaller  modulation  in  the 
channel  conductance. 

§5.  Conclusions 

The  velocity  modulation  field  e fleet  transistor  has  been 
studied  by  self-consistent  particle-field  Monte  Carlo 
methods. 

This  study  shows  that  current  switching  can  be 
achieved  by  the  velocity  modulation  concept. 

The  perpendicular  transport  of  electrons  between  the 
two  channels  is  achieved  in  extraordinarily  short  times 
(  =0.2  picoseconds).  This  is  due  to  short  distances  involv¬ 
ed  between  the  channels  (1000  A),  and  very  high  average 
electron  velocities  (4-6*  101  cm/s)  perpendicular  to  the 
interface.  However,  the  necessary  interface  charge 
readjustment  takes  much  longer  than  perpendicular 
transport,  and  is  accomplished  only  in  times  comparable 
to  the  longitudinal  (source-to-drain)  switch-on  time. 
With  the  proper  choice  of  gate  biases  the  switching  speed 
of  the  VMT's  will  be  superior  to  that  of  the  conventional 
GaAs  field  effect  transistors.  We  estimated  speed  up  of 
abut  1.5  picoseconds  for  1  micrometer  (0.4  micrometer 
gate  length)  device. 
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